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ABSTRACT 

The fraction of Compton thick sources is one of the main uncertainties left in under- 
standing the AGN population. The Swift Burst Alert Telescope (BAT) all-sky survey, 
for the first time gives us an unbiased sample of AGN for all but the most heavily ab- 
sorbed sources A'^h> 10^^ cm~^). Still, the BAT spectra (14 - 195 keV) are time-averaged 
over months of observations and therefore hard to compare with softer spectra from the 
Swift XRT or other missions. This makes it difficult to distinguish between Compton- 
thin and Compton-thick models. With Suzaku, we have obtained simultaneous hard 
(> 15 keV) and soft (0.3 - 10 keV) X-ray spectra for 5 Compton-thick candidate sources. 
We report on the spectra and a comparison with the BAT and earlier XMM observa- 
tions. Based on both flux variability and spectral shape, we conclude that these hidden 
sources are not Compton-thick. We also report on a possible correlation between excess 
variance and Swift BAT luminosity from the 16 d binned light curves, which holds true 
for a sample of both absorbed (4 sources), unabsorbed (8 sources), and Compton thick 
(Circinus) AGN, but is weak in the 64 day binned BAT light curves. 

Subject headings: X-rays: galaxies, galaxies:active 



1. Introduction 

The Burst Alert Telescope (BAT) is a very hard X-ray (14 - 195 keV) telescope surveying the 
entire sky in search of Gamma Ray Bursts. Between bursts, it conducts an all-sky survey. Since the 
BAT is sensitive abo ve 14 keV, it is un biased towards the absorption which greatly affects surveys 



in other wavebands (jMushotzkvl |2004| ) . Analysis of the properties of an unbiased active galactic 
nuclei (AGN) sample will allow us to discover heavily absorbed sources with little or no optical 
indication of AGN emission. This is important because it will give us an indication of the true ratio 
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of absorbed/unabsorbed AGN in the local universe. Since the fraction of C ompton thick so urces 
is one of the main uncertainties left in understanding the AGN population (jComastril |2004| ) . the 
fraction of Compton thick BAT sources will indicate, at the very least, limits on the flux and spectral 
shape for local Compton thic k populations. Also, s ince AGN are li kely the major contributors to 
the cosmic X-ray background (IComastri et al.lll995l : lGilli et al.lll999l ). the distribution of properties 
of the unbiased BAT AGN sources will allow for more accurate models of the X-ray background. 

This work concentrates on five AGN sources detected in the 9-month BAT AGN survey, which 
were associated with 4 sources in the BAT catalog (one pair of sources, NGC 6921 and MCG +04- 
48-002, is associated with SWIFT J2028. 5+2543), w ith a detection sign ificance of A.8-a or above. 
There are 153 AGN sources in this 9- month catalog (jTueller et al.ll2008l ). all with BAT fluxes (14- 
195 keV) in excess of 10~^^ ergs~^ cm~^ and an average redshift of 0.03. In this paper, we present 



Suzaku follow-ups for interesting sources from a previous XMM-Newton follow-up study. 



In 



Winter et al.l (|2008l ). we presented the X-ray properties from ~ 10 ks XMM-Newton obser- 
vations of 22 previously unobserved (in the X-rays) AGNs from the BAT 9-month catalog. Of 
these sources, half were absorbed sources with Nn> 10^3 cm . We found that 5 sources had 
spectra consistent with Compton-thick sources (iVH> 1.5 x lO'^^cm"-^) and 4 more were "hidden" 
AGN, sources wit h lar ge scattering f r action s and a low ratio of -Fo.5-2fcey/-^2-iofcey as discussed in 
Ueda et al.l (j2007l ) and lWinter et al.l ([20091). However, with the 10 ks observations these heavily 
absorbed sources did not have enough counts to distinguish whether a reflection model or a power 
law model better fit the data. Further, adding the BAT data to better constrain the fit is po- 
tentially problematic since the BAT spectra are time averaged over months of observations. Only 
with Suzaku observations could we obtain a simultaneous hard (> 15keV) and soft (< lOkeV) 
spectrum, necessary for further investigation of the complex absorbed sources. 

This paper is the second in a s eries of papers pres enting Suzaku follow-ups of Swift BAT- 
detected AGNs (following paper one: lEguchi et al.l ([2003)). In this paper, we present results of our 
analysis of the Suzaku XIS and HXD PIN spectra for 5 heavily obscured sources (NGC 1142, Mrk 
417, ESO 506-G027, NGC 6921, and MCG +04-48-002). In § [21 we present the data analysis. The 
spectral analysis, focusing on variability (§ [3|) and interesting features (§ d]), follows . In § [H we 



summ arize our results. Throughout our analysis, we use the wilm ISM abundances ([Wilms et al 
2OO0l ). 



2. Data Analysis 



Based on successful proposals during Suzaku AO-1 and AO-2, we obtained observations of 
NGC 1142, Mrk 417, ESO 506-G027, NGC 6921, and MCG +04-48-002 in the HXD nominal 
pointing mode. All of these so urces were absorbed {uh > 10^^ cm~^) , with complex spectra in our 
previous XMM-Newton study (jwinter et aklbood ). In Iwinter et al.l \20o4 ). we noted that MCG 
+04-48-002 and NGC 6921 combined are likely the BAT source, since both of these AGN are within 
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the 6'error circle of BAT. In our XMM-Newton analysis, NGC 6921 was the brighter source with 
a pn count rate of 0.255 ctss^^ compared to 0.096 ctss~^ for MCG +04-48-002. However, in the 
Suzaku observation, NGC 6921 is extremely dim with an XISl count rate of 0.003 cts s~^ compared 
to 0.015 cts s-i for MCG +04-48-002. With a field of view of 1°, the HXD PIN spectrum of this 
observation is also likely a combination of both AGN. However, given the faint state of NGC 6921 in 
the observation, this spectrum is probably dominated by emission from MCG +04-48-002. Details 
of the Suzaku observations for NGC 1142, Mrk 417, ESQ 506-G027, NGC 6921, and MCG +04- 
48-002 are shown in Table [H An analysis of the AO-1 observation of NGC 1142 is also included in 



Eguchi et al.l (|2008l ) (p aper 1 in this series). For details of the analysis of the XMM-Newton and 



Swift BAT spectra, see IWinter et al.1 ((20081). 

To extract each of the Suzaku spectra, we used the cleaned version 2.0 processed event files 
supplied by the Suzaku team. For processing of the XIS data, we combined the 3x3 and 5x5 edit 
modes for the front-illuminated (FI), XISO and XIS3, and back-illuminated (BI), XISl, CCDs in 
XSELECT. The source spectra were extracted from a circular 20" region, centered on the source. 
Background spectra were extracted from 40" regions located nearby regions free of emission. Re- 
sponse matrices (rmf) and ancillary response matrices (arf) were then generated using the newest 
versions of xissimrmf gen and xissimarf gen. Following this, we combined the two XIS-FI spectra 
with MATHPHA and the response files with addrmf and addarf . The spectra and rmf files were re- 
binned from 4096 to 1024 channels, which is still larger than the r esolution of th e CCDs, in order to 
slow down the time the standard spectral fitting software, XSPEC (jArnaudlll996l ). takes to load and 
fit the spectra. The combined FI spectra and XISl spectra were each grouped with the response 
files with 20cts/bin, using grppha. 

For the HXD data, we only used data from the PIN instrument. The spectra were ex- 
tracted using XSELECT for both the PIN event file and the corresponding instrumental tuned 
background file, supplied by the Suzaku team. Both spectra were generated using the same 
good time intervals. The source spectrum was further corrected for instrument dead time. In 
addition to the instrumental background spectrum, we generated a cosmic X-ray background in 
XSPEC using the model suggested by the Suzaku team (a cutoff power law of the form CXB{E) = 
9.412 X IQ-^ X {E/lkeVy^-'^^ x exp(-^/40A;eF) photons cm'^ g-i FOV"^ keV"^). We then com- 
bined the instrumental and cosmic X-ray backgrounds with MATHPHA. Finally, the spectra were 
binned, using grppha, with the background and standard response files from the Suzaku CALDB to 
a signal-to-noise ratio of 3 a. 



3. Variability 

Compton-thick sources are AGN where our line-of-sight to the source is blocked by obscuring 
matter that has an optical depth of r > 1 towards Compton scattering (A'^h> 1-5 x 10^^ cm~^). At 
these optical depths, much of the emission is reflected and not direct emission. For Compton-thick 
sources, the column densities are so high that little to no direct emission escapes below lOkeV. As 
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such, one possible marker of a Compton-thick source, whose emission is dominated by reflection, 
is a lack of variability. For instance, analysis of the recent 100 ks Suzaku observation of the mildly 
Compton thick source NGC 4945 (A^h~ 5 x 10^^ cm~^) found no variability below 10 keV, consistent 
with a reflection dominated spectrum, while the high energy spectrum (> lO k eV), where the direct 
emission is not obscured, showed a factor of two variability (lltoh et al.ll2008l ). 



Further, one of the closest Compton-thick AGN, located in the nucleus of Ci rcinus, is observe d 
as consistently not variable. C i rcinu s has been o bserved extensive l y with ASCA (IMatt et al.l 19961 ) , 
XMM-Newt on dMolendi et ahlbood ). Chandra (jSambruna et al.l boOll : iMassaro et al.1 bood ). and 
Beppo-Sax (jGuainazzi et al.lll999l ). Between these observations, which span approximately 9 years, 
there is no significant change in the shape of the spectrum or the flux of the source. Similarly, the 
recent Suzaku observation of Circin us also shows a r eflection d ominated spectru m with a flux 
consistent with previous observations (jYang et al.ll2008l ). However, I Yang et al.l ([2003) note possible 
low level (2-15%) variability in the spectra both above and below lOkeV, whose uncertainty is 
high due to uncertainties in the pin background (the source is only 7% above the background) and 
possible contamination from ULX sources. If our target sources are Compton-thick, like the AGN 
in Circinus, we would expect little variability in spectral shape and flux. However , if our sources 
resemble "changing-look" sources, sources s uch as NGC 136 5 ([Risaliti et al.ll2005l ) which change 
for Compton-thick to Compton-thin stages (jMatt et al.ll2003l ). we would expect to find variability 
in the absorbing column. 

There are two types of variability that we can probe with our Suzaku observations. These 
types include variability on the time scale of the observation (short term variability) and variability 
between observations for an individual source. In this section, we first characterize short term 
variability in the Suzaku observations with XIS. With fewer counts in the HXD spectra, we do not 
perform a similar analysis at higher energies. To probe variability between observations for individ- 
ual sources, we compare variability in spectral parameters, including spectral index, flux, column 
density, and the fluorescent Fe-Ka line, between our previous 10 ks XMM-Newton observations and 
the Suzaku XIS observations. Finally, we compare variability between the HXD PIN and BAT for 
NGC 1142, the only source for which we have two HXD PIN observations. 



3.1. Variability during the Observations 

3.1.1. Light Curve Analysis 

To test for variability during the individual Suzaku observations, we first constructed 128 s 
light curves (0.1 - 12keV) for both the source and background regions used to extract the spectra 
using XSELECT. Light curves were computed for XISO, XISl, and XIS3. Computing average count 
rates for each light curve, we find that the XISO and XIS3 source light curves are 4-6 times higher 
than the background rates for NGC 1142, Mrk 417, and ESQ 506-G027. The XISl rates are only 
1.5 - 1.6 times higher than the background. For MCG -1-04-48-002 and NGC 6921, the source count 
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rates are lower, compared to the background rates. MCG +04-48-002 has average XISO and XIS3 
rates of about 2 times the background with an XISl rate of only 1.25 times the background. NGC 
6921 is dimmer, with XISO and XIS3 rates 1.4 - 1.5 times the background and XISl average rates 
only 1.2 times the background. 

For each of the XIS light curves, we subtracted the background rates to create a net light 
curve. The XISO, XISl, and XIS3 net light curves were then combined to create an average net 
light curve. For each of these average net light curves, we calculated the normalized excess variance, 
a measurei nent of the va r iabilit y amplitude in the light curve, and ^ values. Excess variability, as 
defined bv iNandra et al.l (|l997l ). corresponds to: 



1 



Nil 



N 



(1) 



Here, // corresponds to the unweighted mean, N is the number of points in the light curve, Xi is 
the count rate at z, and Oi is the corresponding er ror in the count rate. We computed errors on 
al^g using equatio n 11 fromlVaughan et a.1.1 (120031 ). Since the excess variability is dependent on 
observation length (jLawrence &: Papadakislll993l ). we divided the long NGC 1142 (observation 1) 
observation into two evenly spaced observations, computing cr^ms ^-iid individually for each half 
of the observation. These variability measurements are shown in Table [2j We also constructed light 
curves binned by the orbital time scale (~ 5760 s), with the results of this analysis also shown in 
Table [21 

Light curves binned by 5760 s (the orbital time scale) are shown in Figure [TJ While little 
variability exists on rapid time scales sampled by the 128 s light curves, the 5760 s light curves do 
show variability for most of the sources. The amplitude of this variability, however, is not very 
high, with changes in count rates spanning < 0.1 ctss^^ in the combined XISO + XISl + XIS3 light 
curves, corresponding to changes on the 60% or lower scale from the lowest to highest count rates 
in individual light curves. The reduced values range from 0.31 - 2.77, with the exception of 
NGC 6921, whose high value (11.18) may be the result of poorer statistics from the higher ratio of 
background count rates to source count rates. Clearly, there is no strong variability (factors of 2 or 
higher) present on these timescales. While low luminosity A GN, like our obscured targ et sources, 
are thought to be more variable than high luminosity AG N (IBarr &: Mushotzkvl 119861), AGN are 



more variable on long time scales than short time scales (jBarr &: Mushotzkvlll986l : 
19971 ) and so the lack of short term variability is not surprising. 



Nandra et al 



In addition to the full 0.1-12 keV light curves, we constructed 5760 s binned net light curves 
for three energy bands: 0.1-3 keV (L), 3-7 keV (M), and 7-12 keV (H). The average count rate and 
values for each observation are recorded in Table [3j We chose these energy bands to separate 
the soft emission, the region including the Fe K-a emission features, and the hard band emission. 
From our analysis, we find that variability in these regions is not the same for all of these "hidden" 
sources. In paticular we find that the observations of NGC 1142 show more variability at the 
highest energies, while Mrk 417, MCG +04-48-002, and NGC 6921 show more variability in the 
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soft emission. Further, the region including the Fe K-a emission, which contains the most counts 
for all observation except that of NGC 6921, shows little to no variability on these timescales. 



3.1.2. 



af.^g-L Relation 



In order to test whether the anti-correlation between luminosity and c'^rns ( Nandra et alJll997l ) 



is also seen in our sources, we needed to obtain the X-ray fluxes from spectral fits to the energy 
spectra. To model the XIS spectra, we simultaneously fit the XIS f ront-i lluminated and back- 
illuminated spectra with the same standard model used in lWinter et al.l (120081 ) (0.3-10 keV spectra). 
This model consists of a partially covered power law spectrum with an Fe K -a feature at 6.4 keV 
(implemented in XSPEC as pcfabs*tbabs(pow+zgauss) const). The partial covering model is a 
multiplicative model in which the direct AGN emission is partially blocked by material in the line 
of sight, defined as: 

M{E) = / X e-^H'^(^) + (1 - /). (2) 

Here A'nis the absorbing column in the line of sight (in units of atoms cm~^), (j{E) is the photo- 
electric cross section, and / is the covering fraction (ranging from to 1). A large covering fraction 
(/ ~ 1) indicates that either much of the direct emission is blocked or that a very small fraction 
of the emission is scattered into our line of sight. Galactic ab sorption from the Milky W ay is 
accounted for with a second neutra l abso rber (tbabs), fixed to the lDickey LockmanI (Il990l ) value 
(listed in Table 1 of I Winter et al.l (j2008l ) along with z, Seyfert type, and host galaxy type). The 
Galactic A'h values were obtained using the web version of the A^h f tool. Results of these fits 
are presented in Table [H Due to the low signal-to-noise in the NGC 6921 spectrum, the error bars 
are not constrained and encompass the full range of values. Therefore, w e fixed the power la w 
component to 1.75 (the average value of F obtained for the BAT AGNs in I Winter et al.l (|2009l )) 
and record the upper limit on column density and scattering fraction. 

As Table m shows, the values from this model are not optimal. The residuals to the model 
show features not accounted for by this simplified model, for instance there is a soft component 
evident in the spectra of NGC 1142. To investigate these features, we present our detailed spectral 
analyses in § HI For our current study of variability, this simpler partial covering model is sufficient 
to determine luminosities and is further used to provide a direct comparison to the ~ 10 ks XMM- 
Newton spectra (which do not have simultaneous spectra > 10 keV as in the Suzaku spectra analyzed 
in §11]), investigated in the following subsection (§ 13. 2p . 



Computing 0.5-2 keV and 2-10 keV unabsorbed luminosities, using the redshift values recorded 
in Table [H we compare the luminosities with the excess variance measurements in Figure [2l These 
plots show no correlation between normalized excess variance and luminosity, in either band. How- 
ever, the least luminous observations (NGC 6921 an d MCG +04-4 8 -002) do have the highest a^^^ 
values, consistent with the anti-correlation found bv lNandra et al.l (|l997l ) and others. 
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3.2. Variability Between the XMM-Newton and Suzaku Observations 



In lWinter et alj (|2008l ). we found that 13/16 of the sources with SWIFT XRT observations in 
addition to the XMM-Newton follow-ups varied (in flux, power law index, or column density). In 
that study, we had very few counts in the XRT spectra of heavily absorbed sources (< 60 counts 
for sources with A'^h> lO^^cm"^), making it difficult to compare spectra for the most absorbed 
sources. However, we are now able to determine whether our five AGN target sources are variable 
between the XMM-Newton and Suzaku observations. 

To quantify differences in the spectra, we fit the 0.3-1 keV XIS spectra w ith the same partial 
covering model used for the XMM-Newton observations in lWinter et al.l (j2008l ) 13. ip . The power 
law index for NGC 6921 wa s difficult to constr ain and so we again fixed this value to the average 
AGN photon index of 1.75 (jWinter et al.l 12003 ). In Figure O we plot the column density derived 
from the partial covering model, spectral index (F), and 6.4 keV equivalent width versus observed 
2-10 keV luminosity for all of the observations. These plots reveal a number of results. First, we 
find that there is a great change in 2-10 keV luminosity for both NGC 6921 and MCG +04-48-002. 
The luminosity of MCG +04-48-002 increases by one magnitude between the XMM-Newton and 
Suzaku observations (~ 1 year apart), while NGC 6921's luminosity drops by two magnitudes. 
None of the other three sources vary to such a degree, however, all show signs of variability in 
luminosity between observations. Therefore, the heavily obscured sources are varying on time 
scales of ~ 0.5-1.5 years (the time between the Suzaku and XMM-Newton observations). 

Second, it is clear that Ny{^ F, and Fe K-q EW are higher during the lower luminosity ob- 
servation for individual sources. An anti-correlatio n between Fe K EW and l uminosity is known 
as the X-ray Baldwin/ "Iwasawa-Taniguchi" effect ( Iwasawa &: Taniguchi 1993ll and ha s been seen 



in a number o f AGN samples, for instance in the radio quiet samples of iJiang et al.l (|2006l ) and 
Bianchi et al.l (|2007l ). Thus, we would expect to see this anti-correlation, ev en amid our h eavily 
obscured sources. A correlation has also bee n no ted between F and l uminosity (jPorquet et al.ll2004; : 
Piconcelli et al.ll2005l : IShemmer et al.ll2006l ). In lWinter et al.l (|2009l ). we did not find a correlation 
among the 9-mo nth BAT AGN samp le. However, we noted that this correlation was seen for indi- 
vidual sources in IWinter et al.l ([2003) and therefore suggests that sources do have higher spectral 
indices at higher luminosities. It is unclear why we see the opposite effect in our obscured sources 
(i.e. an anti-correlation - note that NGC 6921 is the exception, however, F was fixed to the average 
value of 1.75 in the lower luminosity observation). It is possible that this is the result, in part, of 
unconstrained error bars. In particular, this could be the case for MCG +04-48-002, whose spectral 
parameters were difficult to constrain in the lowest luminosity observations. This may also be the 
cause of the anti-correlation seen between iVnand luminosity, which is dominated by NGC 6921 and 
MCG +04-48-002. For NGC 6921, however, the column density is poorly constrained. Further, a 
value of iVH= 3 x 10^^ cm^^ is beyond the limitations of the tbabs model, since it does not treat 
multiple scatterings, important for A'^h>> 1 x lO^^cm"^. Another possible cause of the observed 
anti-correlation between F and luminosity could be the lack of our inclusion of a refiection model. 
It is possible that refiection plays an important role for some of these sources, because of this, we 
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will discuss models including reflection in our detailed fits in § [H 

To further test variability between the XMM-Newton and Suzaku observations, we simulta- 
neously fit the XMM-Newton pn and Suzaku front-illuminated spectra. Initially, we fixed all the 
parameters to the same values. Next, we allowed the model parameters (A'h and covering fraction 
from pcf abs, F and its normalization, and a constant value to account for flux differences (const 
in xspec)) to vary, recording the change in x^. For all of the observations, we used the partial 
covered power law model presented in § 13.11 However, for NGC 1142 we also added a model for 
collisionally-ionized diffuse gas (apec with kT = 0.63 keV and a normalization of 0.0021) and an 
Fe edge &t E = 7.07 keV with r = 0.40. Each of these additional models signiflcantly improved 
the fit (Ax^ = 178 for apec and 129 for zedge) to the NGC 1142 spectra. In Figured we plot 
the error contours for the 99%, 90%, and 68% confidence levels between F and A'^nfor the simul- 
taneous Suzaku (2) and XMM-Newton fits to the NGC 1142 spectra. As the figur e shows, both 



iVn a nd F are constrained to a fairly narrow parameter space. Also, as we noted in lWinter et al 



(|2008l ). the simple model of a partially covered power law yields fiat power law indices (F ~ 1) 
which are not necessarily physical. Spectral fits with more complex refiection models were not well- 
constrained with the XMM-Newton spectra in the 0.3-10 keV and thus the motivation for obtaining 
these Suzaku observations was to extend the energy range for spectral fits so that we could test 
more complex and possibly more physical models (i.e. reflection). We present the results of these 
fits, which constrain the power law indices to steeper, more physical values, in § HI 

Results of the simultaneous fits are presented in Table El Accompanying plots of the unfolded 
energy spectra are presented in Figure El In the table, in addition to including the A^^ values for 
allowing the model parameters to vary for the different observations, we include a measurement 
of the observed fiux variability. The statistic used is {Fmax — Fmin) / Favg^ which was computed 
for both the soft (0.5-2 keV) and hard (2-10 keV) band. These measurements show that the most 
variation is seen in MCG +04-48-002 and NGC 6921, in both hard and soft bands, and in the hard 
band for Mrk 417. The variability in MCG +04-48-002 is purely in the observed fiux - the shape 
does not change significantly, as indicated by small A^^ for changes in absorption, the apparent F 
from the partial covering model, and Fe K EW. However, both the shape and flux of NGC 6921 
change drastically, with a very signiflcant Ax^ in A'^h(99) and apparent F (47). Finally, the hard 
flux and apparent F (Ax^ = 249) change appreciably in Mrk 417. Here, the flatter F (~ 1.0) 
corresponds to the higher flux observation (the Suzaku observation). From Figure El it appears 
that the difference in shape is most notable below 2 keV. 

In Figure [6l we plot the ratio of our statistic {Fmax — Fmin)/Favg for the 0.5-2 keV/2-10 keV 
bands versus the measurement in the hard band (2-10 keV). The values comp uted for our obscure d 
target sources (circles) are compared to the unobscured sources (square) from I Winter et aP (j2008l ). 
From this plot, we find that the ratio of soft to hard variability tends to be lower for obscured 
sources. Among our sources, only ESO 506-G027 shows more variability below 2keV than above. 
One explanation for the difference in the ratio of soft band to hard band variability between obscured 
and unobscured sources is that the soft band light for the obscured sources is not primarily direct 
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AGN emission, in agreement with lBianchi et al.l (|2006l ). Additional sources of soft X-ray light (e.g. 
binaries, ionized gas) can dilute the AGN signatures at low energies and thus reduce the amplitude 
of variability. 



In lWinter et al.l ( 20091 ). we discussed the fact that emission from X-ray binaries, star formation, 
or hot ionized gas could cause the soft emission measured for the complex AGN. The 0.5-2 keV 
luminosity of all of o ur sources lies below 2 x 10^^ ergss"^, within the range of emission expected 
from star formation ([Ranalli et al.l l2003l ) . To test this further, we added an apec model to the 
simultaneous spectral fits discussed. The apec model is a model for collisionally ioni zed emission 
which can also mimic the emission from photoionized gas (jGuainazzi &: Bianchil 120071 ) . Therefore, 
the apec model can be used to model soft emission that is non-AGN (starburst), AGN (photo- 
ionized gas), or both. We found that this model can account for the soft emission in all of the 
sources except ESQ 506-G027, for which the reduced value becomes unacceptable at 2.3. Since 
the apec model well describes the spectra of the remaining sources and the soft emission is not 
strongly variable, this suggests that for all but ESO 506-G027 the soft emission is dominated by 
non-AGN emission (i.e. star formation or hot gas) , AGN emission (photo-ioni zed gas), or both. 
Photo-ionized emission is extended (e.g., ISako et al.l (j200ll ): iBianchi et al.l (j2006l )) and expected to 
show no variability on time scales of years or less. Since the observations of NGC 1142 showed no 
variability in the soft emission (from the 5760s binned light curves), it is plausible that the soft 
emission in this source is from photo-ionized gas. 

Another result we find is that the amount of hard variability is lower in the unobscured sources 

((Fmax — Fm.in) I V--- < 



MFg^ig ^. 0.5) than the obscur ed target sources (see Figure [6]). While our results on Sy 



Is (jWinter et al.ll2008l ) agreed with those of iNandra et al.l (|l997l ). finding more soft band (< 2keV) 
variability, the obscured sources are different. As we have shown, the source of this variability is 
both from change in flux (for instance in MGG +04-48-002) and shape (i.e., changing A'^h and/or 
apparent T). Since t he obscured sources in the BAT sample tend to have lower luminosities than the 
unobscured sources (jWinter et al.ll2009l ). it is likely that the greater 2-10 keV va riability is the result 
of the sources having a lower luminosity. Therefore, this result is consistent with lBarr &: Mushotzky 
(|l986l ) who first showed an anti-correlation be tween variability ampli tude and X-ray luminosity. 
This result is also consistent with the results of lBeckmann et al.l (120071 ) who found the 14-195 keV 
BAT light curves of bright obscured AGN to vary more than unobscured bright AGN. 



3.3. Variability Above 10 keV 

3.3.1. Swift BAT light curves 

To determine whether our so urces vary above lO keV, we first examined the publicly available 
400 day Swift BAT light curves (jBaumgartner et al.l 120081 ). We chose to bin the light curves by 
16 days, roughly half a month, to investigate variability on this timescale, which corresponds to 
the shortest time scale for which all of the sources are well detected. Binning the light curves by 
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16 days, we computed the average count rate, excess va r iance , and x value for our target sources in 



addition to the unobscured sources from lWinter et alj (120081) . for comparison. We also include the 



14-195 keV luminosity for each source, from lTueller et al.1 (|2008l ). Figure [7] shows the 16 d binned 



light curves for the 4 Swift BAT AGN sources. Results of our analysis are recorded in Table [6j 
While the reduced values of 0.99, 1.12, and 0.84 indicate little variability in the 14-195 keV band 
for ESQ 506-G027, MCG +04-48-002/NGC 6921, and Mrk 417 (NGC 1142 is much higher at 3.87, 
indicating significant variability) , Figure [7] clearly does indicate some variability in all of our target 
sources. Similarly, the unobscured comparison sources show the same variability as the obscured 
sources (indicated by similar cx^^g and reduced measurements, as well as the light curves which 
are not shown). 

To illustrate the similarity in 0"^^^ between the obscured and unobscured sources, we plot 
excess variance versus 14-195 keV luminosity in Figure El Clearly, the obscured and unobscured 
sources occupy the same regions in both luminosity and excess variability. Interestingly, there is a 
correlation between excess variability and luminosity for these sources (obscured and unobscured). 
Fitting the data with an ordinary least squares bi-sector line, we find: 

log ol^, = (1.70 ± 0.48) X log Li4-i95 + (-75.4 ± 20.8). (3) 

Computing several statistics to determine the significance, we find a coefficient of determination 
i?^ = 0.43, a Spearman rank correlation coefficient = 0.76, and Kendall rank tau correlation 
coefficient of r = 0.58. This shows that the relationship of cr^^j, oc xi-70±o.48 ^j^g 14-195 keV 
band is statistically significant. Further, from an analysis of the BAT light curve of the Compton 
thick source Circinus, with /dof = 0.75/23 for a constant source and an excess variance of 
•^rms — 2-8 X 10~'^, whose BAT luminosity is logL = 42.07, we find that this heavily obscured 
source also follows this relation. While the similar relationship between obscured and unobscured 
sources indicates that the same physical mechanism creates the emission in the 14-195 keV band, it 
is unclear why the specific 0"^^., oc i^i-70±o.48 relationship exists, but more data points are necessary 
to increase the statistical significance of this relation. To test whether this relationship holds with 
different binning sizes, we also computed the excess variance and values for 64 day bins, the 
results of which are also included in Table [H While a similar relationship exists between excess 
variance and 14-195 keV luminosity, with cr^^^ a -^14-195^^) the significance of this relationship is 
weak (with similar i?^ of 0.42 but = 0.26 and Kendall rank r = 0.26). Assuming that the 14- 
195 keV is composed of both direct and reflected emission, one possible explanation for a possible 
correlation between luminosity and variability in the 14-195 keV band is that the more luminous 
sources have a smaller contribution from reflection, which may be constant on these time scales. 
However, since the correlation is weak, particularly in the 64 day light curves, it is uncertain whether 
the correlation would be statistically signiflcant in the entire BAT sample. 
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3.3.2. Comparison of Suzaku pin and Swift BAT spectra 

One of our main goals in obtaining Suzaku observations of our targets was to obtain simulta- 
neous spectra above and below 10 keV. This is especially important si nce time averaged 1 3-montli 



BAT spectra have been created for the 9- month BAT AGN sources ([Tueller et al.ll2008l ). While 
these spectra exist, since they are time averaged we have no way of knowing how the spectrum 
changes with time in the 14-195 keV band. As our analysis of the BAT light curves indicates, the 
sources do vary in flux on the ^ half a month time scales investigated. To investigate how the 
shape and flux change for individual sources, it was imperative to obtain multiple Suzaku pin ob- 
servations for direct comparison with the BAT spectra. Towards this end, we obtained two Suzaku 
observations of NGC 1142 (from our AO-1 and AO-2 observations). 

To test the variability or lack thereof in the > lOkeV spectra of NGC 1142, we simultaneously 
fit the pin spectra and BAT spectrum with a simple absorbed power law model (tbabs*pow). 
With no variation between the spectra, this model yields a poor fit with A^h= 1-52 x lO^^cm"^, 
r = 1.72, and x^/dof= 230/61. The residuals from the fit show that while the first Suzaku 
spectrum and the BAT spectrum are not badly fit, the second Suzaku observation is. The fit can 
be improved by allowing the column densities to vary (x'^ = 76/59 when NniBAT) = 2.8 x 10^^, 
A^H(Suzakul)= 2.1 x 10^^, and A''H(Suzaku2)= 9.17 x lO^^cm"^) or allowing the power law index 
of the second Suzaku observation to flatten (x^ = 82.7/59 when A^h= 6.3 x lO^^cm"^, T (BAT, 
Suzakul)= 1.83, and T (Suzaku2) = 0.66). The best statistical fit, = 67/59, is obtained when 
the fluxes of the spectra are allowed to vary (shown in Figure [9]), with A'h= l.Olj^Q^ x 10^^ cm~^ 
and r = 1.53+°;2^. The measured fluxes in the 15-50 keV band are 4.4 (BAT), 3.3 (Suzakul), and 
1.6 (Suzaku2) xlO^^^ ergs~^ cm~^. 

Thus, we clearly find flux variation from the Suzaku pin spectra of NGC 1142, the only source 
for which we have multiple observations. This variation can be well modeled as flux variability 
alone, not statistically requiring variations in the spectral shape. However, with few data points, 
particularly in the dimmer Suzaku 2 observation, we can not rule out additional variations. It 
is important to note that the analysis of the Swift BAT light curve of NGC 1142 found it to be 
the most variable source (x^ = 3.87). However, from our analysis of the XIS light curves and 
XMM-Newton and XIS energy spectra we found NGC 1142 the least variable source in the 0.3- 
10 keV band. Therefore, this source is not necessarily a typical example of variability in the BAT 
AGN catalog. It does show, though, that it is important to allow for the flux normalization to 
vary for the BAT spectra when jointly fitting the BAT spectra with other X-ray spectra (i.e. from 
XMM-Newton or Suzaku). 



3.4. Discussion of Variability 



As already mentioned, one possible indicator of a Compton thick source is a lack of long-term 
variability. The AGN source in Circinus is one of the closest and well-studied Compton thick 
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sources and the luminosity and spectral shape has not changed appreciably over the decade it has 
been monitored in the X-ray band. This is clearly not the case for our five "hidden" AGN. While 
there is no evidence of strong short term variability during the Suzaku observations (~40 - 80 ks), 
there is clear variability between the XMM-Newton/Suzaku observations which are separated by 
~ 1 - 1.5 years. This variability is seen in both the spectral parameters (i.e. A^H; T (though with 
large error bars on the T from the partial covering model this is not statistically significant)) and 
flux. The flux variability is found both above and below lOkeV, with a change in the NGC 1142 
pin flux by a factor of two and order of magnitude changes in the 0.3-10 keV fluxes of NGC 6921 
and MCG +04-48-002. 

From our analysis of the variability in the Swift BAT 14-195 keV light curves (16 day and 64 
day time scales), we found that there is potentially no difference between unobscured and obscured 
sources, a result that must be quantified with an analysis of the entire Swift BAT AGN sample. If 
this holds for the entire sample, it suggests that the same physical mechanism underlies both types 
of sources. Additionally, an analysis of the light curve of Circinus shows that this Compton thick 
source is also consistent with the correlation we found between variability and BAT luminosity in 
the sample (see Figured]). As we suggested, the possible correlation between 14-195 keV luminosity 
and variability may indicate that the relative contribution from a reflection component is lower at 
high luminosity and higher at low luminosity. Since Circinus has a lower 14-195 keV luminosity and 
is reflection dominated, it supports our hypothesis. However, the correlation is weak, particularly 
in the 64 day binned light curves. 

At softer X-ray flux (below 10 keV), the unobscured and obscured sources do differ. In Figure[3l 
we showed that F tends to be steeper /higher in the lowest luminosity observation of an individual 
source. This contradic t s our results comparing XMM-Newton and Swift XRT spectra of unobscured 



AGN in iTueller et al.l (|2008l ). We also found that unlike unobscured AGN, our obscured sources 
have larger scale variability on average in the hard X-ray band (2-10 keV) than the soft X-ray 
band (0.5-2 keV) (Figure [6|). This illustrates a key difference between the obscured and unobscured 
sources. One possible explanation for this difference in variability in the soft band is that the 
soft emission seen in the obscured sources is not direct emission. Given the connection between 



nuclear star forming regions and Seyfert 2 nuclei (IVeilleuxll200ll ) in addition to the lower soft band 



luminosity, it is possible that the soft emission in the obscured sources is dominated by galactic 
emission. However, the fact that some variation is seen shows that there likely is a component from 
the AGN as well in the soft band. 

From our analysis of NGC 1142, we found that the variability above lOkeV is not correlated 
with the variability below lOkeV. In fact, the flux barely changed in the 2-10 keV band (Figured]) 
while there is a factor of 2 change in the Suzaku pin spectra (Figure [9]) . One possible explanation is 
that there is a time delay between the 14-195 keV and 2-10 keV bands, where the softer X-rays are 
the result of some type of reprocessing (i.e. reflection) of the direct emission seen in the 14-195 keV 
band. If this is the case, observing the source in the accompanying low 2-10 keV flux state would 
place constraints on the physical location of the reprocessing material. However, given the anti- 
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correlation between T and luminosity in the 2-10 keV band (for obscured sources), there may be a 
more complex interplay between direct and reprocessed emission in a clumpy absorbing region. 

4. Detailed Spectral Analysis 

In the previous section, we conducted a detailed study of variability in the Suzaku observations, 
between the Suzaku and XMM-Newton spectra, the Swift BAT light curves, and between the Suzaku 
pin and Swift BAT spectra. This study helped us to conclude that the soft emission (< 2keV) 
in all but ESO 506-G027 is not necessarily direct AGN emission (from the lack of variability and 
good statistical fit with an ionized gas model). We also found that the emission above lOkeV is 
variable and that in order to include the time-averaged BAT spectra in our spectral fits, we must 
allow the flux normalization to vary. The presence of variability in all of our sources suggests that 
these AGN are not heavily Compton thick (r >> 1 towards Compton scattering). 

In this section, we examine the spectral signatures of the sources, in particular looking for the 
contribution of reflection for the sources with Suzaku HXD pin data. Since reflection is the main 
signature of a Compton thick source, this will provide further clues to the nature of these sources. 
In addition to looking for refiection signatures, we provide detailed analysis of the spectra for all 
of our sources, characterizing the high signal-to-noise Suzaku observations of the 5 hidden AGN. 

4.1. Suzaku XIS Spectral Fits for NGC 6921 and MCG -^04-48-002 

With an angular separation less than the 6' resolution of Swift's BAT and 1° field of view of 
Suzaku's HXD, both the BAT and pin spectra are combinations of these two X-ray sources. Thus, 
we can not use these spectra to constrain the high energy portion of their spectra. For NGC 6921, 
the low number of counts (330 counts in XISl) for this source, due to the lower flux in the Suzaku 
observation, can not be used for a more detailed analysis than that presented in § [3l Here we 
found that the column density was not well constrained (log Ah ~ 25, a column density above 
the limitations of tbabs, and with errors encompassing the entire allowed range of Ah). It is 
possible that the column density of the source is heavily Compton thick in this observation, i.e. 
log Ah ^ I/ct preventing transmission of non-reflected emission. However, as in the previous 
XMM-Newton spectrum, a strong Fe K-a EW, a signature of reflected emission, is not present. 
Thus, the spectra of this source continue to be a puzzle. 

For MCG +04-48-002, the superior spectral resolution of the Suzaku XIS detectors allows 
us to better constrain the parameters of the observed Fe K-a feature. Extracting spectra of the 
calibration source, Fe-55, we found that the errors on detectable line widths range from 0.013 keV 
- 0.024 keV. Therefore, we can detect line widths above this range. To determine the line width 
and energy of the Fe K-a line, we fit the combined XIS0+XIS3 and XISl spectra with the partial 
covering model used in Table [H We extended the energy range to 12keV to better constrain 
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the continuum. Within this energy range, we find that the column density and covering fraction 
are the same as in Table HI within the errors. The spectral index is slightly steeper, measured 
as r = 2.36^Q4g. From the residuals of the spectral fit, we found evidence for an Fe XXV K-a 
emission line in the combined XIS0+XIS3 spectrum (not seen in the XISl spectrum). Adding this 
feature improved by 4.5, showing the line to be marginally significant. We could not constrain 
the width, fixing this to 0.01 keV. From our spectral fit, we find an energy of G.GS^^qq^ keV and EW 
of 70+^^ eV for Fe XXV K-a. We find E = 6.38+[]:[j^ keV, a = 0.05+g:|]^ keV, and EW of 184^^^ eV 
for Fe I K-a. The energies of these lines are perfectly consistent with the known optical redshift 
and identification as fluorescent Fe I Ka and He-like Fe. 



4.2. Joint Suzaku XIS, HXD, and Swift BAT Spectral Fits 



For the remaining sources, NGC 1142, Mrk 417, and ESO 506-G027, we performed joint fits 
of the Suzaku XIS, HXD pin, and Swift BAT spectra. The inclusion of the BAT spectra provides 
an extension of the spectrum beyond the HXD pin high energy limit for our sources. Since the 
14-195 keV flux is not constant over the 13 month period used to construct the BAT spectra, we 
included a constant value which was allowed to vary for the BAT spectra. 

As a first step to our fitting process, we used the partial covering model described in § [3] (a 
partially covered power law with an Fe K-a line) to fit the combined XIS + pin + BAT spectra. 
We fixed the pin flux to be 1.12 times the XIS flux (as described in the Suzaku documentation). 
We also allowed the energy and width of the Fe K-a line to vary. As a next step, we replaced the 
power law model with a cut-off power law - a power law model with a high energy cut-off. As in 
§ [31 we included an apec model for the spectra of NGC 1142. We present the important results 
from these spectral fits in Table [71 One result is that adding a high energy cutoff to the power 
law is statistically significant for all of the spectra (Ax^ > 12). The cutoff energies we measure are 
~ 50-80 keV with error bars of ~ 20 keV, similar to or lower than high energy cutoffs measured 
in some AGNs wit h joint XMM-New t on/Integral spectral fits based on a siin ple absorbed cutoff 
power law model (jMolina et al.l 120061 : [Panessa et al.l [20081 : [Molina et al.[ [2008[ ) . In Figure [lOl we 
plot the error contours (using the steppar command) for the cutoff energy versus power law index 
for the highest resolution observation (NGC 1142 obs. 1). As shown, the computed energy cutoff 
depends on the power law slope (i.e. a fiatter slope gives a lower cutoff energy). Other models 
examined below give higher cutoff energies or show no strong constraint on cutoff energy (i.e. the 
reflection model and a double partial covering model). 

Another important result is that none of the measured columns are in the Compton-thick 
regime (A^h> 1-4 x lO^'^cm"^). Adding the high energy data constrains the columns to lower 
values, though still highly obscured (logA'n > 23). Finally, the measured power law slopes for 
our sources are still flat (F ~ 1-0), much fl atter than the aver age slope of 1.75 measured for the 
entire BAT sample in the 0.3-10 keV band ([Winter et al.[ [20091 ). Since the power law model is an 
approximation to the AGN emission, it is likely that the inclusion of a reflection model or a more 
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accurate Compton scattering model is crucial. 

Before including more complicated models of the AGN emission, we investigate the presence of 
emission and/or absorption features around the Fe K-a emission line. To quantify the importance 
of these additional features, we include the Ax^ values on adding an Fe I K-a line as well as 
additional lines often seen in high signal-to-noise spectra (Fe XXV Ka, Fe I K/3, and Ni I Ka) in 
Table[8l Also, the errors on each of the parameters for these features are included. Since the pcf abs 
model places the 7.11 keV Fe edge at a redshift of 0, we substituted this model with zpcfabs, a 
model which has an additional parameter to shift this feature to the redshift of the AGN. 

Clearly, the Fe I Ka line is significant (Ax^ > 50) in all of these sources, with the spectrum 
of Mrk 417 showing the weakest line. The width of this line, < EW >~ 67 eV, is consistent with 
velocities from about 1800 - 3000 km s^^. These results agree with those from Chandra grating 
results, which place the origin of the Fe K-a lin e in a region between the narrow (500 - 700 kms~^) 



and broad (3000 - 10000 km s'^) line regions (jYaqoob &: Padmanabhanll2004l ). Of the remaining 



emission lines we fit the spectra with, Fe I K/3 is detected in the spectra of NGC 1142 (both 
observations) and ESO 506-G027 and Fe XXV Ka is only detected in t he spectrum o f ESQ 506- 



G027. Ni I Ka is not detected in any of the spectra. As described in lYaqoob et al.l (|2007l ) and 
elsewhere, the ratio of Fe I K/3/Fe I Ka can be used to constrain the ionization state of the gas. 
From the Suzaku observations of NGC 1142 and ESO 506-G027, this ratio is ~ 10%, consistent 
with neutral iron. However, the errors are large (ranging from 2-50%), requiring higher quality 
data to better constrain these values. 

In addition to these features, the addition of an Fe edge at 7.11 keV (in addition to the edge 
built into the zpcfabs model) was significant {Ax^ = 7.5) in the long observation of NGC 1142. The 
best-fit parameters for the energy and optical depth of this feature (added with a zedge model) are 
E = T.W^q'^q keV and r = 0.364^*^0073. Since the 7.11 keV edge is a possible signature of refiection, 
this could mean that reflection is significant only in the first observation of NGC 1142. Alternatively, 
the Fe edge is also an indicator of the Fe abundance. Using the optical depth (r) cal culated from 



the ed ge model and the cross section for H like Fe at 7.11 keV {a = 5.304 x 10 '^^ cm'^ : iHenke et al 



(jl993l )). we calculated the column density of iron (A'h= t/(^) and divided this by the hydrogen 



column density obtained from the p artial covering mod el. Assuming an ISM iron abundance of 



2.69 X 10 ^ the hydrogen abundance ( Wilms et al.ll2000l ). this yields abundances of approximately 



3.5 solar for NGC 1142 (observation 1). The lack of a detection in the second NGC 1142 observation 
suggests that reflection may be the best explanation for this source. However, if the slope of the 
power law is set to 1.0, as in the longer observation, the upper limit on the iron edge optical depth 
is 0.363, in line with the value obtained for the flrst observation. 
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Alternative Models 



While the reduced values from the partially covered cutoff power law models indicate good 
fits to the data (x^ ~ 1-0), the flat power law slopes (F ^ 1.0) are unphysical. A n alternative 



model including a transmitted, scattered, and reflected component was applied by lEguchi et al 



([2003) (Model B) to the Suzaku spectra of "hidden" Swift-detected AGN, including the long ob- 



servation of NGC 1142. This model is represented as tbabsG'ai*(ztbabstr.ans*cutof f plj^ans + 
const*cutof f p l.sf^n,/: + ztbabSrp.//*pexraVrp. f ; + Fe lines) in XSPEC. The parameters of the reflection 
model (pexrav (iMagdziarz &: Zdziarskilll995l )) include photon index, folding energy, abundance, in- 
clination angle, reflection scaling factor (R = ^}/2tt, which we allowed to vary between -5 and 
(negative value s in this mod e l acco unt for the reflected spectrum alone)), and the normalization. 
In Model B of lEguchi et al.l (|2008l ). the normalization of the reflected component is fixed to the 
same value as the transmitted and scattered component (whose spectral indices are also fixed to 
the same value, along with the F value f or the reflected corn ponent). We fix the folding energy to 
300 keV (the default value, also used in lEguchi et al.l (j2008l )). abundances to solar, and cosine of 
the inclination angle to the default of 0.45. In Tabled we include the best fit paramet ers using this 



reflect ion model. While the first observation of NGC 1142 was alre ady included in lEguchi et al 



(1200a), we in clude our own fit, which u ses the wilm ISM abundances (lEguchi et al.l (120081) use angr 



abundances (lAnders Grevessdll989l )) and is binned differently than in 
bin the spectra by 50 cts/bin. For all of the sources, our fits include the Fe I Ka line. We note that 
using the reflection model, we can not constrain a cut off energy when we allow the cutoff energy 
to be a free parameter. We found that the model prefers no cutoff energy (the value floated to the 
upper limits imposed), showing that there is a strong trade off of reflection versus cutoff energy. 
Further, whether the cutoff energy is set at a lower value of 100 keV or 300 keV, we find that the 
photon index from these fits does not change appreciably and is much steeper than the values of 
F ~ 1.0 obtained with the partial covering model. 

The results of this reflection model (Table [9]) indicate good fits to the spectra with reduced 
X^ ~ 1.0 and more "normal" power law slopes (< F >= 1.76). These po wer law index values are 
more typical of the slopes found for the BAT sample (jWinter et al.ll2009l ). However, the reflection 
component measurements, R, show the sources to have mild reflection components at most (with 
the possible exception of ESO 506-G027). A value for \R\ > 1 , as found for the long observation 
of NGC 1142, is unphysic al, as described in lEguchi et al.l (|2008l ). Our value is similar to that from 
the analysis presented in lEguchi et al.l ([200a), where they classify NGC 1142 as mildly or hardly 
Compton thick and explain the high value of R as t he result of the dir ect emission being completely 
blocked non-uniform material in the line of sight (jUeda et al.ll2007l ). These results suggest that 
none of these "hidden" AGN are Compton thick in the sense that they are not reflection dominated 
and have column densities below 1.5 x 10^^ cm^^ (with the possible exception of ESO 506-G027). 
The classification as not heavily Compton-thick is supported by the fact that the sources do vary 
between observations in flux as well as spectral shape. Particularly, this argument applies to NGC 
1142, whose pin spectrum is nearly three times dimmer in the second Suzaku observation. 
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Another model which is a good fit to the data and does not require reflection is the double 
partial covering model. This model assumes that two partially covering absorbers, with different 
column densities and covering fractions, are within the line of sight to the direct emission. Thus, 
this model includes one unabsorbed cor nponent plus three absorbed components with different 
column densities, contrary to Model C in lEguchi et al.l (|2008l ) which includes two different column 
densities for the transmitted component. In Figure [TH the double partial covering model is plotted, 
which is represented in XSPEC as tbabsGa/*zpcf absi*zpcf abs2*(cutof f pi + Fe lines). Fits to 
the individual spectra are shown in Figure [l2] and the key parameters from this fit are listed in 
Table [TDJ Since this model constrains the cutoff energy to a lower value in the long observation of 
NGC 1142, Ecutoff = 14:3~^j'^]ieV , than in the reflection model where it is unconstrained, we fix 
Ecutoff = 100 keV for the observations where it could not easily be constrained with the double 
partial covering model. 

With this double partial covering mode l, the power law in dices are more consistent with the 
average value for AGNs in the BAT sample (jWinter et al.ll2009l ). Further, with this model we can 
explain the difference between the spectra of NGC 1142 in two ways - a difference in absorbing 
columns as well as a st e eper spectrum i n the secon d, dimmer observation . Contrary to the results of 
Piconcelli et al.l (j2005l ') . JPorquet et al.l (120041 '). and lShemmer et al.l (|2006l ). as well as our own results 
from Seyfert 1 sources (jWinter et al.ll2008l ). the accretion rate or flux is anti-correlated with spectral 
index for the two observ ations of NGC 1142. The s a me is true of the Fe Kg equivalent w idth, 
contrary to the results of llwasawa Sz Taniguchil (|l993l ). iPage et al.l (|2004l ). Ijiang et al.l (|2006l ). and 
others. This highlights another key difference between Sy Is and the hidden sources, suggesting 
th at there likely i s a co mplex environment of varying absorbers. This result is consistent with those 
of iRisaliti et al.l (120021 ) . who concluded that variations in column density are ubiquitous in Sy 2s 
and the product of clumpy absorbers. 



The Difference Spectrum for NGC 114-2 

With two Suzaku spectra of NGC 1142, we created a difference spectrum - a spectrum con- 
structed from subtracting the lower flux observation from the higher flux observation. To do this, 
we used the Ftool mathpha to subtract the observation 1 combined xis 0+3 source spectrum from 
the observation 2 source spectrum. We did the same for the background spectra. We then used 
the Ftools addrmf and addarf to add the response files from the two observations, weighting them 
by their respective exposure times (observation 1 was weighted as 0.688 and observation 2 was 
weighted as 0.312). In the same way, we created the pin difference spectrum between the two 
observations. We binned the spectra in the same manner as in § [2j 



Unlike the difference spectra of other sources, such as MCG-5-23-16 (jReeves et al.ll2007l ). the 
difference spectrum of NGC 1142 is not a pure power law. A partial covered cutoff power law 
is not a good statistical fit to the data, with x^/dof = 953/436. As shown in Figure [T3| there 
is emission seen below 1 keV. This soft emission shows that the soft spectrum did in fact change 
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between the two observations. Adding an apec model to the base partial covered cutoff power law 
model (zpcf abs*cutof f pi) improves the fit by Ax^ = 44. The temperature of this component 
(kT = 0.63]to 04 keV) is consistent with the values we previously found for both observation 1 and 
2. 

In addition to a soft component, emission and absorption features are present in the difference 
spectrum. The Fe I Ka line is very significant, with Ax^ = 195. Both the energy and width of 
this line are consistent with the previous measurements. The equivalent width is the same as in 
observation 1, with EW = 225^32 eV. The Fe I K/3 line is also present, but much less significant 
with Ax^ = 4. Since these lines are seen as emission features in the difference spectrum, this shows 
that the flux of these lines scales with the continuum flux (i.e. the line fiux is higher when the 
continuum is higher). 

Further, three other features are also present. At low energy, an emission line is present at 
E = 1.09 lb 0.02 keV, with Ax^ = 15 and EW = 56 eV. The energy of this line is consistent with 
Fe XXII L (E = 1.053 keV). Additionally, a 7.2 keV absorption feature is significant in the difference 
spectrum (Ax^ = 9). This suggests that the edge energy is changing between observations or that 
there is an absorption feature at this energy which is reacting to changes in the spectrum. The final 
feature seen is another "absorption" feature, with A^^ = 15 and E = 7.67±0.01 keV {EW = 98 eV). 

Upon adding these features, the best fit model to the difference spectrum yields an acceptable 
fit of x^/dof = 670/422 (1.59). The best-fit parameters for absorption and power law index are: 
iVH= 7.7 ± 0.8 X 10^3 cm-^ Cvr = 0.995 ± 0.001, and F = imtoH- Compared to the best-fit 
parameters for observations 1 and 2 (Table [7]), the difference spectrum has a higher column density, 
slightly higher covering fraction, and steeper power law index (particularly in comparison with 
observation 1). 



4.3. Discussion of Spectral Properties 

From our detailed spectral analysis, we found that a reflection component is very weak in the 
spectra of our target sources. Also, none of the sources have a 1 keV EW Fe I Ka line. Since our 
sources do not meet these two criteria for Compton-thick spectral classification, we classify them 
as Compton-thin. All of the spectra are well-fit by a partial covering or double partial covering 
model, implying that multiple absorbing columns lie along our sight line to the hidden sources. 

Through our spectral fits, we find that the derived column densities and power law components 
change depending upon the model. For instance, very fiat power law components are obtained using 
a power law model for NGC 1141, Mrk 417, and ESO 506-G027 (Table [7]) while steeper values are 
found using the reflection model or double partial covering model. This shows that a flat power 
law is not a good indicator of a Compton thick or reflection dominated spectrum and that other 
criteria like a 1 keV EW Fe I Ka emissio n line or lack of lo ng-term variability (except in the case of 



'changing-look' Compton-thick sources (|Matt et al.l 120031 )) should be met. Additionally, we found 
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that the cutoff energy for the power law component is also model dependent. With a simple partial 
covering model, we obtained cutoff energies of f» 100 keV, while the cutoff could not be constrained 
with the reflection model (and was therefore set to 300 keV). Using the double partial covering 
model, we constrained the cutoff energy to ~ 150 keV in the highest signal-to-noise observation of 
NGC 1142. Therefore, caution must be taken in constraining cutoff energies, since these values are 
highly model dependent. 

In Figures [T ^ 115 ^ and 1161 we plot the relationship between N}^, T, and Fe I Kq EW and lumi- 
nosity (absorption corrected in the 2-10 keV and 10-50 keV bands) and Eddington ratio (L/LEdd)- 
The model parameters are obtained from the double partial covering model (Table [TOj) and from 
Table H] for MCG -1-04-48-002 (shown only in the top left panel of each figure) . The Eddington lu- 



minosity was obtained from the 2MASS derived masses shown in lWinter et al.l ( 20091 ) (1.3 x 10 



l38 



{Mbh /Mq) ergss"^). As shown in Figure [T4l there is no relationship between Nn and either lumi- 
nosity or Eddington ratio. However, as shown in Figure \T5\ while there is no apparent correlation 
between F and luminosity, there is possibly an anti-correlation between F and L/LEdd (in both the 
2-10 and 10-50 keV bands). In the variability section, we already mentioned this b ehavior, which 



i s opp o site the correlat i on se e n between F an d lumi nosit y in other samples ( i.e. iPorquet et al 



(1200J); iPiconcelli et al.l (l2005l ): IShemmer et al.l tOO^ ). In IWinter et al.l (l2009l ). we had found no 



correlation of any kind between F and luminosity or Eddington rate but interpreted this as a result 
of not having asample of similar objects (like the PG quasars used in the samples mentioned). In 
Winter et al.l (|2008l ). we did see the correlation for multiple observations of individual unobscured 
sources (comparing Swift XRT and XMM-Newton spectra). If there is an anti-correlation for the 
absorbed hidden sources, this suggests either that the model being used is wrong (even the double 
partial covering model may be too much of a simplification) and/or that there is a different physical 
mechanism responsible for the spectra of these sources. However, since the error bars are very large 
on the F measurements, the anti-correlation is not statistically significant. 

In Figure [TBI we plot the Fe I Ka EW versus luminosity and Eddington ratio. As shown, 
there is no correlation seen between the EW and luminosity. However, there is a slight correlation 
between EW and the Eddington ratios. Fitting a line to each plot, we find: 



log EW = (-0.38 ± 0.16) X logiL^'^f^y LEdd) + (l-U ± 0.50) 



log EW 



-0.32 ± 0.12) X log{L^^^%/LEdd) + (1-38 ± 0.34). 



(4) 
(5) 



While this correlation is not as stati stically signifi c ant a s the correlation we found between EW 
binned by Lg °Tn/ivKrfrf in F i gure 9 of I Winter et aD ^20od ). with R"^ = 0.27 and 0.34, the result is 



consistent. In 



tent with the results of 



Winter et al.1 (l2009ll. w e found EW oc (L^°!To/^£;dd)~°■^^^°■°^ ^hich is also consis- 



Bianchi et al.l (|2007l ) for radio quiet AGN {EW oc (Lboi/LEdd) 



-0.19±0.05\ 



Therefore, the accretion rate and Fe I Ka emission are tied together and the mechanism which 
creates this line is the same between the unabsorbed and absorbed sources. 



- 20 - 



5. Summary 



In this paper, we present an analysis of Suzaku observations of five absorbed sources from the 
9-month Swift BAT AGN catalog, NGC 1142, Mrk 417, ESQ 506-G027, NGC 6921, and MCG 
+04-48-002. These nearby sources {< z >= 0.023) were only recently detected in the X-rays by 
Swift's BAT. Through our study of these sources, we conducted an analysis of both their variability 
properties and detailed spectral properties with Suzaku. 

From a study of variability during the 40 ks Suzaku observations (as well as the 80 ks AO- 
1 observation of NGC 1142), we found little short term variability (i.e. only small amplitude 
variability of < 0.1 cts"^ in the combined XIS0+XIS1+XIS3 hght curves) on 128 s or 5760 s time 
scales. On longer time scales of 0.5-1.5 yrs, between the XMM-Newton and Suzaku observations 
of these sources, variability is present. Both the shape (i.e. through changing A^h and/or T) and 
flux change, most notably for NGC 6921 and MCG +04-48-002 whose fluxe s change by an orde r 
of magnitude or more. Further, this result is consistent with the results of iRisaliti et al.l (|2002l ). 
suggesting that variability is ubiquitous in absorbed sources. Comparing the chan ge in flux of the 
absorbed sources with a sample of unobscured sources from I Winter et al.l (120081 ). we found that 
the obscured sources show significantly higher levels of hard band (2-10 keV) variability. This 
highlights a potentially important difference between absorbed and unabsorbed sources and is in 
agreement with previous findings wh ich show that the soft component in the absorbed sources is 
at a larger scale (jBianchi et al.ll2006l ). 



When we examined variability in the 14-195 keV band, through 400 d Swift BAT light curves, 
we found that both the obscured and unobscured sources showed similar levels of variability. Inter- 
estingly, the 16 d binned light curves showed a correlation between excess variance and luminosity: 
-^14-19^5 which is also found in the 64 day binned hght curves {a'^rns ^ -^14-195 ^^)' 
where the correlation is weaker. This is, however, the first time such a correlation has been noted. 
One clear implication of this analysis is that AGN do vary above 10 keV. This variability does not 
appear to be correlated with variability in the softer bands, which is evident from NGC 1142, the 
least variable source below 10 keV and the most variable above 10 keV. 

From a detailed analysis of the Suzaku spectra of these "hidden" sources, we present a few 
interesting results. Among these, we constrain the properties of the Fe I Ka feature (central energy, 
width, intensity, and equivalent width), finding the width of the line consistent with velocities from 
about 1800 - 3000 km s-^ Thus, our Suzaku results agree with Chandra grating results, placing 
the origin of the Fe K-a lin e in a region between the narro w (500 - 700kms-i) and broad (3000 - 
10000 km s"i) line regions (jYaqoob &: PadmanabhanI |2004| ) . We also detect additional features in 
the spectra, including a significant 7.11 keV Fe edge in the spectra of NGC 1142, Fe XXV Ka in 
the spectra of ESO 506-G027 and MCG +04-48-002, and Fe I K/3 in the spectra of NGC 1142 and 
ESO 506-G027. 



Finally, based on variations in flux, a small reflection component, and Fe I Ka EW << 1 keV, 
it appears that none of our sources are true Compton-thick objects. A more likely explanation for 
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their properties is heavy absorption (logA'n > 23) in a clumpy environment (jRisahti et al.ll2002l ). 
In support of this, we find good fits to the spectra with a double partial covering model, a model 
where two absorbers partially cover the central emission, with no reflection component required. 



L.W. would like to thank Alex Markowitz (UCSD) for advice on generating light curves and 
processing the Suzaku data and Yuxuan Yang (University of Illinois) for access to the Suzaku 
spectrum of Circinus. We also thank Satoshi Eguchi for discussions on the Suzaku data sets. 
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Fig. 1. — Light curves binned by 5760s (orbital period) for NGC 1142 observation 1 (top left), 
NGC 1142 observation 2 (top right), Mrk 417 (middle left), ESQ 506-G027 (middle right), MCG 
+04-48-02 (bottom left) and NGC 6921 (bottom right). The light curves are background subtracted 
and added XISO + XISl + XIS3 curves. 
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Fig. 2. — Normalized excess variance, as a measure of variability, versus unabsorbed 0.5-2 keV (left) 
and unabsorbed 2-10 keV (right) luminosity. The individual sources are symbolized as NGC 6921 
(upside down triangle), MCG +04-48-002 (square), ESQ 506-G027 (circle), Mrk 417 (triangle), and 
NGC 1142 (diamond). There is no statistically significant relationship between a^^g and luminosity, 
however, the least luminous sources have the highest normalized excess variance measurements. 
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Fig. 3. — Comparison of the XMM-Newton (| Winter et al.ll2008l ) and Suzaku spectral fits 13. ip in 
the 0.3-10 keV band. We plot A^h, T, and the Fe K-a EW at 6.4 keV versus L2-10. The individual 
sources are symbolized as NGC 6921 (upside down triangle), MCG +04-48-002 (square), ESO 506- 
G027 (circle), Mrk 417 (triangle), and NGC 1142 (diamond). Dashed lines are used to clearly 
distinguish observations of the same source. From the plots, it appears that A^h, F, and Fe K 
EW are higher at lower luminosities for individual sources. Note that the exception to this in the 
F-L2-10 plot is NGC 6921, where F could not be constrained and was fixed to 1.75 in the lower 
luminosity observation. For this source, Nn is an upper limit and the Fe Ka EW could not be 
determined for the Suzaku observation. 
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Fig. 4. — Shown are the error contours for each of the NGC 1142 spectra simultaneously fit, includ- 
ing the two Suzaku observations and the XMM-Newton observation. The errors were computed 
for the photon index (F) and column density (A^^nin units of 10^^ cm^^). The contour levels shown 
represent 99%, 90%, and 68% confidence levels. 
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Fig. 5. — Unfolded spectra showing the XMM-Newton and Suzaku XIS 0-3 simultaneous fits for 
our target sources, rebinned by signal-to-noise for illustrative purposes. NGC 6921 underwent the 
most change between observations (1 year). The source is 2 orders of magnitude less luminous in the 
Suzaku observation and more highly absorbed. Where the error bars are very large (i.e. extending 
beyond the plotted range), a horizontal line is used to represent the location of the Suzaku data 
point. Similarly, error bars extending below the plotted range in the y-axis are excluded for clarity. 
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Fig. 6. — The ratio of {Fmax — Fmin) / Favg in the 0.5-2 keV band/2-10 keV band versu s the value in 



the 2- 10 keV band for our obscured targets (circle) and the unobscured sources from lWinter et al 



([20081) (square). The line represents values where the soft flux variability and hard flux variability 
measurements are the same. This figure shows that there is more hard band variability in the 
obscured sources. 
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Fig. 7. — The Swift BAT light curves for our targets spanning 400 days and binned by 16 days. The 
hght curves represent NGC 1142 (top left), Mrk 417 (top right), ESQ 506-G027 (bottom left), and 
SWIFT J2028.5+2543 (which includes both NGC 6921 and MCG +04-48-002) (bottom right). All 
of these sources show variability during the BAT observations. The line in each plot represents the 
average count rate for that source. 
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Fig. 8. — Plotted are the excess variability measurements in the 16 d (left) and 64 d (right) binned 
BAT lightcurves versus 14-195 keV luminosity. We plot bot h the t arget obscured sources (circles) 
and the comparison unobscured sources from IWinter et al.l (|2008l ) (squares) . There is no visible 
difference between the obscured and unobscured sources. However, our data does show a correlation 



between excess variability and luminosity in this band of cr^^ns ^ 

L 



L 



2.12±0.45 
14-195 



1.70±0.48 
^14-195 



(16 d) or cjI 



(X 



(64 d). A Spearman rank correlation coefficient of 0.76 indicates that the correlation is 
significant in the 16 d binned light curves, but not in the 64 d binned light curves (where = 0.26). 
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Fig. 9. — Comparison of the Swift BAT spectrum (green, time-averaged over 13 months) with the 
two Suzaku pin observations of NGC 1142 (black = observation 1, blue = observation 2) in the 
15-200 keV band. The power law component (F = 1.53) and column density (A^h= 1-0 x 10^^ cm~^) 
are the same for all of the observations, while the flux is allowed to vary by a constant factor. From 
these unfolded spectra, it is clear that the flux is much lower in the second Suzaku observation 
while the first observation has a spectrum consistent with the time-averaged BAT spectrum. 
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Fig. 10. — Using the partial covering model for the long observation of NGC 1142 (Table [7]), we 
present the 99%, 90%, and 68% error contours for the high energy cuttoff and photon index (F). 
This shows that the cutoff energy is dependent of the power law slope, using this partial covering 
model. A flatter slope yields a lower cutoff energy. 
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Fig. 11. — Plotted is the double partial covering model used to fit the combined Suzaku XIS, pin, 
and Swift BAT spectra, implemented as tbabsG'ai*pcf abs*pcf abs*(pow + Gaussian lines)*const 
in XSPEC. 
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Fig. 12. — Plotted is the observed Suzaku + Swift BAT spectrum fit with the double partial 
covering model, along with the ratio of the data to the model, for each of the observations for NGC 
1142, Mrk 417, and ESQ 506-G027. The double partial covering model (Table lOD, shown, yielded 
similar /dof ~ 1.0 values as the reflection model fits (Table [9]). The Suzaku XIS spectra are 
shown rebin to a signal-to-noise of 10 (NGC 1142 obs-1) or 5 (for the remaining sources). 
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Fig. 13. — Plotted is the difference XIS 0+3 spectrum and difference pin spectrum for NGC 1142, 
created from subtracting the low flux Suzaku observation from the high flux Suzaku observation. 
The model used is a partial covering cutoff power law model (zpcf abs*cutof fpl), with a constant 
parameter to allow for flux differences between the XIS and pin spectra. As seen, this model is not 
an adequate fit with significant residuals at soft energies and a prominent Fe I Ka feature. This 
shows that both the soft emission and the iron line change between observations. 




ymate (ergs 




10"* 10" 



z 



10*' 



Lio-sn (ergs s^^ cm-^) 




Fig. 14. — Plotted is the line of sight column density (A^h) versus the unabsorbed luminosity in 
the 2-10 keV band (top left), the 10-50 keV band (top right), the ratio of 2-10 keV luminosity to 
the Eddington luminosity, and the ratio of 10-50 keV luminosity to Eddington luminosity. There 
is no correlation seen between these values. 
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Fig. 15. — Plotted is the spectral index (F) versus the unabsorbed luminosity in the 2-10 keV 
band (top left), the 10-50 keV band (top right), the ratio of 2-10 keV luminosity to the Eddington 
luminosity, and the ratio of 10-50 kcV luminosity to Eddington luminosity. There is evidence of a 
possible anti-correlation between F and the computed Eddington ratios, however, the errors on F 
are too large for this to be statistically significant. 
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Fig. 16. — Plotted is the Fe I Ka equivalent width versus the unabsorbed luminosity in the 2-10 keV 
band (top left), the 10-50 keV band (top right), the ratio of 2-10 keV luminosity to the Eddington 
luminosity, and the ratio of 10-50 keV luminosity to Eddington luminosity. There is no obvious cor- 
relation in the luminosity plots, but a slight correlation {B? = 0.27 and 0.34, respectively) emerges 
in the Eddington ratio plots. The correl ations agree with th e correlation we found in our binned 
EW versus corrected Lz-io/^sdd plot in I Winter etaP (|2009l ^. with EW oc {L^Tiq / LEdd)'^-^^^^'^^ 
and EW « (L^r'5o/ii^<i<i)"°-''^°-''- 
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Table 1. Details for the Suzaku Observations 



Source 


RA (h m s) 


Dec (°"') 


z 


Obs ID 


Date 


Exp (s)* 


Ct Rate* 


NGC 1142 - 1 


02 55 12.2 


00 11 01.0 


0.028847 


701013010 


2007-01-23 


101630, 80580 


0.060, 0.081 


NGC 1142 - 2 








702079010 


2007-07-21 


40570, 36540 


0.041, 0.020 


Mrk 417 


10 49 30.9 


22 57 51.8 


0.032756 


702078010 


2007-05-18 


41507, 13865 


0.04, 0.022 


ESQ 506-G027 


12 38 54.6 


-27 18 28.1 


0.025024 


702080010 


2007-08-02 


35721, 41761 


0.035, 0.027 


NGC 6921 


20 28 28.9 


25 43 24.3 


0.014287 


702081010 


2007-04-18 


41299 


0.008 


MCG +04-48-002 


20 28 35.0 


25 44 00.0 


0.013900 


702081010 


2007-04-18 


41299 


0.026 



*Exposure time and count rate for XISl and PIN. NGC 6921 and MCG -f04-48-002 are in the same observation, 
however, NGC 6921 is very dim. 



Table 2. Variability in Suzaku Observations in 128s and 5760s Bins 



Source 


<Ct Rato*" 


ct2 b 




xVdoF 


<Ct Rato'' 


a2 b 

^ rms 








128 s Bins 








5760 s Bins 




NGC 1142 - 1.1 


6.93 


27 ±3.5 X 10- 


-6 


1.63/398 


6.94 


6.1 ± 1.4 X lO"'^ 


2.77/10 


NGC 1142 - 1.2 


6.99 


10 ± 1.9 X 10- 


-6 


1.35/398 


6.64 


0.02 ± 1.1 X 10-* 


0.99/10 


NGC 1142 -2 


4.99 


38±1.7x 10- 


-5 


1.77/313 


5.06 


0.5 ±4.2 X 10-* 


1.12/8 


Mrk 417 


5.08 


28 ±1.1 X 10- 


-5 


2.07/324 


4.89 


5.6 ±2.7 X 10-^ 


2.15/10 


ESQ 506-G027 


3.58 


8.5 ± 1.0 X 10 


— 5 


1.44/327 


3.32 


5.0 ±4.9 X 10-* 


0.31/6 


MCG-f04-48-002 


4.61 


140±3 X 10- 


4 


2.45/326 


4.24 


13 ±9.1 X 10-6 


2.11/9 


NGC 6921 


1.35 


220 ±0.008 




2.79/326 


0.89 


23 ±4.3 X 10-2 


11.18/9 



^Average background subtracted count rate (10 ^x ctss ^) for the 0.1 - 12keV band (averaged between the 
XISO, XISl, and XIS3 observations). 

''Corresponding excess variabihty measurements (xlO-'^), as defined in lNandra et al.l (jl997| ). 
'^Reduced value divided by the number of bins (dof). 
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Table 3. Variability in Suzaku Observations in 5760s Bins at Specified Energies 



Source 


<Ct Rate>i,'' 




<Ct Rate>M'* 


xVdofM'' 


<Ct Rate>H* 


xVdofe'^ 


NGC 1142 - 1 


1.14 


1.05/20 


4.20 


0.63/20 


1.34 


2.29/20 


NGC 1142 -2 


1.29 


1.14/8 


2.68 


1.17/8 


0.86 


1.60/8 


Mrk 417 


0.73 


2.74/10 


3.04 


2.18/10 


1.08 


1.12/10 


ESQ 506-G027 


0.56 


0.52/6 


1.99 


1.06/6 


0.66 


0.93/6 


MCG+04-48-002 


1.21 


3.86/9 


2.04 


1.88/9 


0.87 


1.00/9 


NGC 6921 


1.05 


9.18/9 


0.05 


1.24/9 


0.11 


3.37/9 



Average background subtracted count rate (10~^x ctss~^) for the 0.1 - 3keV band (L), 3 - 7keV band 
(M), and 7 - 12keV band (H), all averaged between the XISO, XISl, and XIS3 observations. 
''Reduced value divided by the number of bins (dof). 



Table 4. Spectral Fits (0.3-10 keV) with the Partial Covering Model 



Source Ny; Cvr.'- T Fc EW^ y-Zdof F,s, F^/'' L,s., L 



NGC 1142-1 


O.OZ_g 20 


n qqv+o.ooi 

0.»»< _o.o02 


9 00+0.14 
^•^^-0.06 


233111 


1212.2/797 


0.12, 3.60 


43.78, 43.61 


NGC 1142 - 2 


o rQ+0.88 

o-Oy_0.41 




^■^^-0.26 


305+^? 


257.3/180 


0.10, 2.19 


43.60, 43.42 


Mrk 417 




0.974t-.ry 


^■-■^^-0.02 


126^30 


331.7/314 


0.04, 2.99 


42.55, 43.21 


ESQ 506-G027 




0.9761°;^ 




5101^8 


304.3/236 


0.05, 2.22 


42.44, 43.06 


NGC 6921 


307* 


0.554* 


1.75t 




110.4/83 


0.02, 0.05 


41.26, 39.69 


MCG +04-48-002 




n QQn+0 005 


-1 71+0.31 




244.7/408 


0.08, 2.29 


42.48, 42.80 



^Absorption (in addition to the Galactic value) is modeled with the partial covering model. N^f is in units 
of xlO^^cm^^ while Cvr. is the partial covering fraction. 

^Fe K EW in eV, assuming an Fe K-a line with width 0.01 keV at an energy of 6.4 keV. 

^The units for quoted fluxes are xlO~^^ergs~^cm~^. Fs is the observed 0.5-2.0 keV flux and Fh is the 
observed 2.0 - 10.0 keV flux. 

^Logarithm of the absorption corrected (unabsorbed) luminosity in the 0.5-2.0 keV (Lg) and 2.0-10.0 keV 
{Lh) bands. 

* These parameters are upper limits. 

'^This parameter was fixed to the indicated value. 



Table 5. Results of Simultaneous XMM-Newton and Suzaku Spectral Fits (0.3-10 keV) with the 

Partial Covering Model 



Source 


Soft Var.i 


Hard Var.^ 




r2 


Fe EW^ 


NGC 1142 


0.37 


0.52 


47.1 


38.4 


2.9 


Mrk 417 


0.55 


0.80 


23 


249.4 


0.3 


ESQ 506-G027 


0.65 


0.53 


77.6 


-12.5 


22.5 


NGC 6921 


0.88 


1.98 


98.84 


46.7 


0.1 


MCG +04-48-002 


1.10 


1.60 


3 


12.2 






^Observed {Fmax — Fmin) / Favg in the soft band (0.5-2keV) and hard 
band (2-10 keV). 

^Ax^ when absorption (A^jjand covering fraction for the pcf abs model), 
r and its normalization, or the Fe K normalization (zgauss) are allowed 
to vary between the observations. 



Table 6. Variability in the BAT lightcurves (16 d, 64 d) 



Source 


<Ct Rate>'' 


2 b 
^rms 


xVdof 


<Ct Rate>^ 


2 b 
^rms 




-L14-195 






16 d 






64 d 






NGC 1142 


1.46 


67 + 7.3 X 10"* 


3.87/23 


1.42 


120 + 1.3 X 10"^ 


11.51/5 


44.17 


Mrk 417 


0.51 


210 + 0.83 


0.84/23 


0.41 


81+4.5 X 10"=* 


0.39/5 


43.95 


ESQ 506-G027 


1.28 


100 + 0.01 


0.99/23 


1.15 


6.3 + 3.6 X 10-* 


1.49/10 


44.28 


MCG -h04-48-002/NGC 6921 


1.14 


24 + 2.0 X 10"^ 


1.12/23 


1.10 


10 + 1.0 X lO"'^ 


1.02/5 


43.45 


Mrk 352 


0.50 


38 + 0.30 


1.49/23 


0.45 


38 + 5.8 X 10"^ 


1.74/15 


43.27 


ESQ 548-G081 


0.68 


1.6 + 7.7 X 10"^ 


1.14/23 


0.76 


7.6 + 8.2 X 10-^ 


0.82/5 


43.19 


ESQ 490-G026 


0.65 


300 + 0.09 


1.38/23 


0.69 


26 + 3.5 X 10"^ 


1.45/5 


43.71 


2MASX J09043699+5536025 


0.29 


420 + 0.43 


0.87/23 


0.30 


110 + 4.3 X 10"^ 


0.28/10 


44.03 


MCG +04-22-042 


0.52 


4800 + 4000 


1.18/23 


0.57 


4.5 + 1.1 X 10"^ 


1.12/10 


43.99 


UGC 06728 


0.43 


9.3 + 0.19 


1.17/23 


0.53 


0.02 + 3.8 X 10"^ 


1.09/15 


42.54 


WKK 1263 


0.41 


60 + 59 


1.22/23 


0.40 


260 + 0.58 


0.32/5 


43.58 


MCG +09-21-096 


0.58 


58 + 9.8 X 10"^ 


1.37/23 


0.58 


26 + 1.0 X 10"=* 


1.91/5 


43.72 



''Average count rate (10~*x ctss~^) for the 14 - 195 keV band BAT li ghtcurves. 
''Corresponding excess variability measurements (xlO~^), as defined in iNandra et al.l ( 19971 ). 
'^Reduced value divided by the number of bins (dof ) . 

Logarithm of the 14-195 keV luminosity, from lTueller et al.l (|2008h 



Table 7. Spectral Fits (0.3-195 keV) with the Partial Covered Cut-off Power Law Model 



Source Nh'' Cvr.^ T Ecutoff (keV) Ax^*" xVdof 



NGC 1142 - 1 


4 

^••J^-o.io 


0.970; 


f0.004 
-0.003 


n 70+0.20 
' •J-0.11 


32.8+303° 


35.4 


914.0/830 


NGC 1142-2 


4 61+0-29 


0.975; 


fO.006 
-0.007 


-■^•+"-0.12 


80.6+«2?-3 


12.5 


470.1/313 


Mrk 417 


o oQ+0.34 


0.981; 


f0.003 
-0.011 


-■-•^^-0.17 




23.3 


361.5/336 


ESO .")0(i-G027 


"•-'-'^-(liil 


0.979; 


f O.OIl") 
-(1. 010 


1.22+0-:l! 


o7.4t-:;: 


2<S.O 


.S24.0/2G8 



Absorption (in addition to the Galactic value) is modeled with the partial covering model. 
Nh is in units of x lO^^ cm"^ while Cvr. is the partial covering fraction. 

^Ax^ between the partially covered power law and partially covered cut-off power law 
model. 



Table 8. Emission Features with the Partial Covered Cut-off Power Law Model 



Component 


Parameter^ 


NGC 1142 - 1 


NGC 1142 - 2 


Mrk 417 


ESO 506-G027 


Fe I Ka 


Ax' 
Ec 
cr 
/ 

EW 


531.9 

fi 004+0-008 

0.054««1^ 
2471^^ 


173.7 
6.394tnS? 
0.064irol^ 

511? 
367tlt 


54.3 
6 qs6+" "^^ 

O.OOD_o.036 

201^ 
162lfo 


231.2 

ono+0.023 

0.057+ 

79+10 

'^-18 

45ll??i 


Fe XXV Ka 


Ax' 

Ec 

I 

EW 


0.7 

6.722 
0.5+ 
17t 




6.67* 
6t 
29+ 


0.5 
6.67* 

3+ 
22+ 


6.8 

6.6571°:!°^ 
8.21^:^ 


Fe I K/3 


Ax' 

Ec 

I 

EW 


5.0 
7 noq+o osfi 

311^ 

list 


5.3 

R QQQ+0 055 
"•^^^-0.061 

4.911? 
45l?I 


0.1 
7.06* 
13+ 
146+ 


4.4 
7.0481°;-^ 
5.711? 

qQ+32.9 
^^-31.2 


Ni I Ka 


Ax' 

Ec 

I 

EW 


1.7 
7.451 
8t 
63+ 


0.6 
7.45* 
6+ 
66+ 


1.3 
7.45* 
7+ 
124+ 


0.4 

7.45* 
11+ 
99+ 



^The parameters from fitting tfie feature with a gaussian model. These parameters include the Ax^ 
value from adding this component, the central energy of the component in keV, the intensity in units of 
10~^ photons cm~^ s~^, the width of the line (cr, where not indicated this is fixed to O.OlkeV), and the 
equivalent width (eV) with respect to the absorbed cutoff power law model. 

*This value is fixed. 

+ This value is an upper limit. 



- 44 - 



Table 9. Spectral Fits (0.3-195 keV) with the Reflection Model 



Parameter^ 


NGC 1142 1 


NGC 1142 - 2 


Mrk 417 


ESO 506-G027 


r 

fscat (xlO-2) 

N^f^f^ (xl023cm-2) 
R 

kTapec{keV) 

punabs ( X ;^o-" ergs s^i cm-2) 


ii.lO_;^ 43 
1 09+009 

9 cc+o.ee 
—9 on+1-83 

^•^"-0.85 

71+0 04 

' -^-0.06 

1.99 


7.98t°:^t 

1.91^°;?? 

U.OO_Q Q5 

4.55t 
61+° °^ 

"■"-^-0.06 

1.24 


5 55+"-''2 

1 4t^+0-13 
r. 1 7+0.06 
c 94+2.23 

-0.19t 
1.06 


8.4i;?-p7 

-1 C9+0.10 
n 10+003 

1.4lt 
-0 46+0-33 

1.70 


-F'lo-to ( X 10"" ergs s"i cm-^) 


4.48 


1.54 


2.31 


3.00 


xVdof 


898.5/882 


349.7/321 


322.0/342 


263.8/242 



^The parameters from a reflection model represented as tbabsG'ai*(ztbabSfraris*cutof f plf,.a„s + 
const*cutof f plscat + ztbabsre/i*pexravre/; + Fe lines) in XSPEC. The transmitted, scattered, and re- 
flected power law components are fixed to have the same spectral index (F) and cutoff energy (300 keV) . 

^Upper limit on the indicated parameter. 



*Note that in Eguchi et al.l (j2008l ) fscat is the constant factor recorded here divided by the reflection 
factor, R, for values of i? > 1.0. 
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Table 10. Spectral Fits (0.3-195 keV) with the Double Partial Covering Model 



Parameter** 


NGC 1142 - 1 


NGC 1142 - 2 


Mrk 417 


ESO 506-G027 


Nu^ (xl023cm-2) 




7.961°;?^ 


q 7-1+0.38 
"5-' ■'--0.76 


1 9 q-r +1-28 


Cvr^ 




0.9781°:°°! 




0.969t°;°}^ 


A^H^ (xl023cm-2) 
Cvr^ 


q 77+1.59 
' ' -1.36 


0.88tJ:°^ 


^•^-'--0.33 


1.44l°:f5 




0.68li°:°f7 


642+°-^^^ 


0.816tn^^ 


r 


■'■•*J'J_o.l8 


1 D9+0.15 


1 44+0-14 


i-68iS:}^ 


Ecutoff (keV) 


143173 


100* 


100* 


100* 


kTapgQ (keV) 


U.00_o.o3 


62+°-°^ 








"•^"-0.06 


0.12t 






punabs (xl0-"ergss-icm-2) 


2.35 


1.23 


1.05 


1.90 


Fi™-lo ( X 10-" ergs s'^ cm-^) 


4.66 


1.63 


2.32 


3.10 


xVdof 


898.5/881 


338.9/319 


316.1/338 


242.7/240 



^Absorption (in addition to tlie Galactic value) is modeled with the partial covering model. Njj is in units 
of xlO^^cm"^ while Cvr. is the partial covering fraction. Our double partial covering model is represented 
in XSPEC as: tbabs*(apec + zpcf abs*zpcf abs*zedge*(cutof f pi + Fe hues)). The photo-ionized model 
(apec) and Fe edge at 7.11 keV are only used where they are statistically significant. 
Upper limit on the indicated parameter. 



*The value could not be constrained and is fixed to 100 keV. 



